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Abstract 

An effect generated by the nonexponential behavior of the survival 
amplitude of an unstable state at the long time region is considered. 
It is known that this amplitude tends to zero as t goes to the infinity 
more slowly than any exponential function of t. Using methods of 
asymptotic analysis we find the asymptotic form of this amplitude in 
the long time region in a general model independent case. We find 
that the long time behavior of this amplitude affects the form of the 
instantaneous energy of unstable states: This energy should be much 
smaller for suitably long times t than the energy of this state for t of 
the order of the lifetime of the considered unstable state. 
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1 Introduction 

Unstable states \u) of the system under considerations are characterized by 
their decay law, V u (t), 

V u {t) = \a{t)\\ (1) 
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where 

a(t) = (u\u{t)), (2) 

is the probability amplitude of finding the system at the time t in the initial 
state \u) G T~t prepared at time t = 0, Tt is the Hilbert space of states of 
the system, || \u) || = 1 so a(0) = 1, and \u(t)) G Ti. solves the Schodinger 
equation 

ih^\u(t))=H\u(t)), \u{0)) = \u), (3) 

where H denotes the total selfadjoint Hamiltonian for the system. From 
basic principles of quantum theory it is known that the nondecay amplitude 
a(t), and thus the decay law V u {t) of the unstable state \u) decaying in the 
vacuum, are completely determined by the density of the energy distribution 
uj(s) for the system in this state pQ, 

r+°° i , 

a(t) = / u(e) e~n £t de. (4) 



where u>(e) > 0. 

Note that (j4j) and (J2]) together with the condition that \u) is a normalized 
vector mean that there must be 

/+oo 
u(e)de=l. (5) 
-oo 

From the last property one concludes that u{e) is an absolutely integrable 
function, u(e) G L\{— oo, oo). So the amplitude a(t) is the Fourier transform 
of u}(e) (see (j4])) and thus from the Riemann-Lebesgue Lemma it follows that 
a(t) must tend to zero as t — > oo [HE]. 

A condition, which is necessary on physical grounds, that H has a spec- 
trum bounded from below, Spec.(H) = [E min , +oo), and E min > — oo, re- 
duces the set of functions lo(e) G L\{— oo,oo) to such u(e) G L\{— oo, oo) 
that u>[e) = for e < E min and co(e) > for e > E min . Thus in fact the 
integration in (jlj), ([3]) and in similar formulae is taken over e G [E min , +oo). 

In (3] assuming that the spectrum of H must be bounded from below 
and using the Paley-Wiener Theorem [1] it was proved that in the case of 
unstable states there must be 

\a(t)\ > Ae~ bt \ (6) 
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for \t\ — > oo. Here A > 0, 5 > and < q < 1. This means that the decay 
law "P u (£) of unstable states decaying in the vacuum, flTJ, can not be described 
by an exponential function of time t if time t is suitably long, t — > oo, and 
that for these lengths of time V u (t) tends to zero as t —>■ oo more slowly 
than any exponential function of t. The analysis of the models of the decay 

■Y°t 

— n 

processes shows that V u (t) ~ e ft , (where 7" is the decay rate of the state 
\u)), to a very high accuracy for a wide time range t: From t suitably greater 
than some T ~ t = but T > to C^m(£) has nonexponential power-like 
form for short times £ G (£o>^o) ~~ see ; e -g- [31 El El) up to f > r u = -| 
and smaller than £ = t as , where t as denotes the time t for which the long 
time nonexponential deviations of a(t) begin to dominate (see eg., [3], [5] 
|10j). From this analysis it follows that in the general case the decay law 
V u (t) takes the inverse power-like form t~ x , (where A > 0), for suitably large 
t > t as S> t u [3], [5] - [?]. This effect is in agreement with the general result 
([6]). Effects of this type are sometimes called the "Khalfin effect" (see eg. 

mi). 

The problem how to detect possible deviations from the exponential form 
of V u (t) in the long time region has been attracting attention of physicists 
since the first theoretical predictions of such an effect [12J, [13l E] • The tests 
that have been performed over many years to examine the form of the decay 
laws for t ^> t u have not indicated any deviations from the exponential 
form of V u (t) in the long time region. Nevertheless, conditions leading to 
the nonexponetial behavior of the amplitude a(t) at long times were studied 
theoretically [2] - [22] . Conclusions following from these studies were applied 
successfully in experiment described in [23], where the experimental evidence 
of deviations from the exponential decay law at long times was reported. This 
result gives rise to another problem which now becomes important: if and 
how the long time deviations from the exponential decay law depend on 
the model considered (that is, on the form of ui(e)), and if (and how) these 
deviations affect the energy of the unstable state and its decay rate in the 
long time region. 
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2 General long time properties of the nonde- 
cay amplitude 

Many fundamental and general model independent properties of the nonde- 
cay amplitude a(t) follow from the fact that a(t) is the Fourier transform of 
an absolutely integrable function uj(e). So, if one assumes that all derivatives 
u/ fc )(e), (k = 0, 1, . . . , n), exist and that u^(e) G L\(— oo, oo) for all these 
k, (where u/°)(e) = uj{e)), then analyzing general properties of the Fourier 
transforms it is easy to find that for t — > oo 

Ht)\ < p (7) 

where < C < oo (see [2]). So if the derivative of u>(e) exists and LO^\e) G 
L\{— oo,oo) then taking into account relation (jHJ) the following estimation 
follows 

j > \a(t)\ >Ae~ htq (8) 

as t — > oo. 

Much more information about asymptotic properties of a(t) being the 
Fourier transform of u(e) one can find if the assumption required by physics 
that the spectrum of H is bounded from below is used, eg. by e = E min > 
— oo. This means that from this time we will consider only such u)(e) G 
L\{— oo, oo) that u(e) = for e < E min and u(s) > for e > E min . In this 
case there is 

a(t) = / u{e) e h et de. (9) 

J E 

instead of the formula (jlj) for the amplitude a(t). Using methods of the 
asymptotic analysis it is not difficult to find an asymptotic expansion for 
large values of t of the Fourier integral of this type [2H [25l [26] . 

Let us consider for a start relatively simple case when \im £ ^ Emiri+ u(e) = f 
cu > 0. Let derivatives Lo^ k \e), (k — 0, 1, 2, . . . , n), be continuous in [E min , oo). 
(that is let for e > E min all u^ k \e) be continuous and all the limits 
lim e ^£; m . n+ exist) and let all these be absolutely integrable 

functions then, as can be easily verified, (see [24" t [25l [26] ). 



fc=0 
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where u { = lirn £ _> jBmin+ a/ fc >(e). 

Let us now consider a more complicated form of the density to(e). Namely 
let oj(e) be of the form 

u(e) = (e - E min ) x 77(e) G Li(-oo,oo), (11) 

where < A < 1 and it is assumed that r)( k >(e), (k = 0, 1, 2, . . . , n), exist 
and they are continuous in [E min ,oo), and limits lim e _»_E min+ rj( k \e) exist, 
lim^oo (e — E min ) x rj^ k \e) = for all above mentioned k, then 



a{t) 



ih . —i vp. . + 

\ e h^minl* 

t 



a. 



:(t) + ( 



— ) Oin-xit) 



1 ifr\2 , s 

+ (- —) a n _ 2 (t) 

/ ih,s , , 

+ (" ~) On- 3 (*) + ... 



(12) 



where (compare [25l [26] ) 

n ^r([+A) -iiLtet?) 

z=o 



and % J = \im E ^ Emin+ ^(e), r]^(e) = r](e) and j = 0, 1, . . . , n. 

The main difference between asymptotic expansions ( fTUl) and ( fl~2l) is that 
the amplitude a(t) in (JTUl) is obtained as the Fourier transform (J9j) of such 
uj(e) that co>(e) = for e < E min and co(E min ) > whereas the expansion 
f|T2|) is the asymptotic expansion of the Fourier transform (Q for another 
type u;(e): namely for u(e) such that uj(E min ) = (see (TTT]) ) . 

From ffTUl) . ffl~2l) it follows that the estimation ([8]) is true for all physically 
admissible 00(e). This estimation holds for t — > 00, strictly speaking for 
£ > t as . For t ~ r u , t < t as the nondecay amplitude a(t) takes the form 

a(t) ~ e _i i( E °~t^), (14) 

to a very high accuracy [3] , [5] [11] , [H] - [22] . In this formula E® denotes 
the measured energy of the unstable particle described by the state-vector 
\u). There is El > E min . 
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For the energy densities 00(e) leading to the asymptotic form of the am- 
plitude a(t) of type (TTOT) the time t as can be found by comparing the square 
of the modulus of the amplitude a(t) from the relation f|T^|) and the square 
of the modulus of the leading component of ffTUl) . So t as can be found by 
solving the following transcendental equation 

e~ 1 f t = tf( u{E ™ n) ) 2 . (15) 

This means that the value of t as depends on the model considered: it depends 
on the density uj(e) and on the 7°. If oj(e) has the form (TTTT) then a similar 
method can be used to find a corresponding equation for t as . 



3 Energy of unstable states at long time re- 
gion 

Using the Khalfin's estimation ([6]) of the decay law ([T|) one can examine the 
asymptotic properties of the decay rate j u of an unstable state \u). In a 
general case the decay rate 7^ equals, 

,.x def h &P u (t) 

7 « = 7 " (t) = - vjt) -&T- (16) 

From (|T6|) . ([!]) and (|T4|) one infers that 

lu = lu {t)= 1 l (fort < to.), (17) 

which is obvious. From ([6]) it follows that in the asymptotic case t — > 00 the 
decay rate 7«(t) can not be larger than [27] . 



7„(t) ~ bgr^, (18) 

t— »oo 

where /i = 1 — g > 0. So for £ > t as one finds that 7«(^)|^ 00 < 7° for every 
physically admissible u>(e), and in general that lim t ^oo j u (t) = 0. 

The problem is how the energy E u of the unstable state \u) behaves for 
t > t as . The solution of this problem follows from the observation that the 
amplitude a(t) can be found either by solving the Schrodinger equation (jSJ) 
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or using the equation for the projection of the state vector ( see [27J and ref- 
erences one can find therein), which in the case of one-dimensional subspace 
7i\\ of states 7i spanned by the normalized vector \u) has the following simple 
form 

ih ^_ = K a ( t ) j a(0) = 1; (19) 

where h u is the "effective Hamiltonian" for the one-dimensional subspace 
of states Ti\\. In general, h u can depend on time t, h u = h u {t) [27J 128] - 
One meets this effective Hamiltonian when one starts with the Schrodinger 
Equation fl3]) for the total state space TC and looks for the rigorous evolution 
equation for the distinguished subspace of states TC\\ C TC. There are many 
approximate methods to calculate h u [27] but taking into account the problem 
raised above their use is not necessary. It is sufficient to use the property 
that the exact effective Hamiltonian h u (t) must fulfill the following identity 

K = h u (t) ^ift^-L. (20) 
at a{t) 

Direct application (1201 to the relation (114j) yields 

h u (t) = h° u = E u - l -~f° u , (fort<U, (21) 



which could be expected. Note that from (11611 and (1201) it follows that simply 

lu (t) = -2Q(h u (t)). (22) 

Similarly, the real part of h u {t) is the instantaneous energy, E u (t), of the 
system in the state \u) under considerations 

E u = E u (t) = M(h u (t)). (23) 

(Here 3? (z) and 3 (z) denote the real and imaginary parts of z respectively). 

Note that relations (TSOj) and f[T9l establish a direct connection between 
the amplitude a(t) for the state \u) and the exact effective Hamiltonian h u (t) 
governing the time evolution in the one-dimensional subspace TCu 3 \u). 
Thus the use of the relation ( |20|) is one of the most effective tools for the 
accurate analysis of the early- as well as the long-time properties of the 
instantaneous energy and decay rate for a given qausistationary state \u(t)). 
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Now let us analyze the asymptotic properties of h u (t) for t — > oo. For 
the densities u(e) leading to asymptotic expansion (flU]) one finds that for 
t — > oo, 

era £ 



X 



o (!) q ( 2 ) z'^ 2 

(3) /*^x3 



-4^(y) +•••)• (24) 

The next step is to use the relation ([2D]) . So one should now divide ff24l) by 
( TTUj) and then collect together all components of the same order with respect 
to ( j ). As the result one obtains the asymptotic form of h u (t) for t — > 00, 

One obtains a similar form of h^(t) for the amplitude a(t) given by formulae 
( |T2l) . ( |T3l) . Indeed starting from ( fl2l) one finds after some algebra that 

Kit) = E mm -c^-c 2 (^) 2 - c 3 (|) 3 + . . . , (26) 

where Ci, 02,03,... are complex numbers with negative or positive real and 
imaginary parts. 

A surprising conclusion following from the result (125]) is that in the long 
time region the leading component of the asymptotic form of the decay rate 
j u (t) has the same form for a large class of physically admissible models, 

lu(t) = TuWUoo = "29 (W)) * 2 5 . (27) 

Using the relation (125]) we find that in the long time region the instantaneous 
energy E u (t) takes the following form 

Eu(t) d = f E u {t)\ t ^ = R(h™(t)) 

- E mm - ^(j) 2 + ... ^ i& (28) 
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As one can see for all densities to{e) such that u> = uj{E min ) > the long 
time properties of the leading components of the energy E u (t), contrary to 
the properties of the decay rate 7«(t), depend on the density co(e). The result 
( 128]) seems to be even much more surprising than ff27|) . 
We have 

lim E?(t) = E min . (29) 

t— >oo 

Note that the same result follows from (126]) which means that the relation 
( 129]) is a model independent. Taking into account that E® > E min the 
following conclusion follows: for every model (that is for every U)(e) ) there 
exists such > t as that 

E?(t) < El (for t > t^). (30) 

Note that results (j2"51) - ([3"U]) are purely quantum effects and that they follow 
from basic assumptions of quantum theory. 

4 Final remarks 

The problem if the long time deviations from the exponential form of the 
decay law affect the energy of the decaying state has been studied in [2] using 
a model defined by oj(s) having a form of the truncated Lorentz function and 
assuming that E min = 0. It is easy to verify that inserting into (fTUl) . (|25|) 
E m in = and u(e) used in [9] reproduces relations obtained there. Results 
obtained in Sec. 2 and Sec. 3 show that long time behavior of the decay law 
V u (t) as well as the effective Hamiltonian h u (t) do not depend on a specific 
form of the density uj(s) but they depend rather on general integral and 
analytic properties of the density u(e). From these results it follows that for 
all u(s) having the same integral and analytic properties the amplitudes a(t) 
and the effective Hamiltonians h u (t) have the same long time behavior. So 
these results generalize and complete essentially analysis performed in [S]. 

The estimation (jSj) of \a(t)\ for t — > oo follows from basic assumption 
of quantum theory. Similarly, estimations (flQl) . ([121) of a(t) at long time 
region are obtained using only very general assumptions on the form of the 
energy density u(e) and the fundamental assumption that there exists a 
minimal energy E min > — oo in the system under considerations (i.e. that 
Spec.(H) = [E min , +oo)). So, they should hold for every physical system 
fulfilling these general requirements. In general, as it follows from the analysis 
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performed in Sec. 2, two types of the long time asymptotic expansion of the 
amplitude a(t) can be observed depending on the continuity properties of the 
density u(e) at the point e = E min . If the density lo(e) is a discontinuous 
function of e at e = E min : uj(e < E min ) = and u(E min ) > 0, then the long 
time asymptotic of a(t) is given by the formula ffTUl) . A particular, typical 
example of such u(e) is the truncated Lorentzian distribution function lul^e), 

Me) = =-e(e- E mm ) J" (31) 

2vr ( e _ Elf + (^) 2 

where AMs a normalization constant, > E min , and 0(e) = {1 for e > 
0, and for e < 0}. This distribution is the basis of many studies 
of decaying systems (see, eg. [21 ED])- Many unstable systems will have an 
initial state energy distribution u(e) that is close to u)l(s) for all e values. 
It has been proved [5] that in such cases the decay law V u (t) for the system 
and decay law V^(t) resulting from ujl{s) must be close to each other for 
all values of t. In general the density ou(e) having Lorentz (Breit-Wigner) 
shape is known from the response of a harmonically bound elektron with 
a dissipative term, models of resonance behavior and many other physical 
problems. On the other hand, if lo{s) has the form ffTTj) . that is if it is 
continuous at e = E m i n , then the long time form of the amplitude a(t) is given 
by the relation f|T2|) . A particular case of this type density distribution is 
oo(e) which can be found when one considers short-range potential models of 
quasi-stationary states: One can find such a density for finite-width barriers 
as well as delta barriers and with or without a potential inside the barier, 
etc. (see, eg., [151 ESI 13 I2Q1 [211 122] and references one can find therein). In 
general in the models mentioned densities uj(e) are proportional near E m i n 
to the square root of the energy e, 

uj{e) ~ y/(e-E min ), (32) 

for e > E min , (i.e. they correspond with A = | in (fTTj) and (|T2|) ), and usually 
there is E min = in these models. The another example of the density of 
this type is the density io{e) obtained when one considers the decay of an 
unstable particle into two particles [7]. The form of long time asymptotic 
expansions for a(t) can differ from expansions ffTUj) and f|T2|) for u(e) being 
discontinuous at a point (or some points j — 1, 2, . . . ,) e — Ej > E min . 

The source of the effect described by the relation (1251) is the long time be- 
havior of the amplitude a(t). The relation (12"U|) establishes a direct connection 
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between the properties of the amplitude a(t) and the properties of the instan- 
taneous energy E u {t) and decay rate j u (t) of the unstable \u) at the instant 
t, (see f l2Tj) . fl23|) . fl22|) ). A possibility observe the "loss of energy" described 
by relations ff28l) . fl30l may arise while trying to test the long time properties 
of the nondecay amplitude a(t) after a suitable modification of such tests. 
So, considering a possibility of a suitable modification of the test described 
in [23] in such a way that the emitted energy (frequency) of the luminescence 
decays could be measured which could make it possible to test relations (|28|) . 
( pOj) seems worthwile. In general, all these long time properties of unstable 
states should not be expected to have an effect on laboratory processes but 
it seems that they can affect some long time astrophysical processes. 
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